Descriptive Thermodynamics 
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Thermodynamics (in concert with its sister discipline, statistical physics) can be regarded as 
a data reduction scheme based on partitioning a total system into a subsystem and a bath that 
weakly interact with each other. Whereas conventionally, the systems investigated require this form 
of data reduction in order to facilitate prediction, a different problem also occurs, in the context of 
communication networks, markets, etc. Such "empirically accessible" systems typically overwhelm 
observers with the sort of information that in the case of (say) a gas is effectively unobtainable. 
What is required for such complex interacting systems is not prediction (this may be impossible 
when humans besides the observer are responsible for the interactions) but rather, description as a 
route to understanding. Still, the need for a thermodynamical data reduction scheme remains. In 
this paper, we show how an empirical temperature can be computed for finite, empirically accessible 
systems, and further outline how this construction allows the age-old science of thermodynamics to 
be fruitfully applied to them. 
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MOTIVATION: THE PASSIVE DEMON 

In the following we outline a program for an alterna- 
tive approach to thermodynamics and statistical physics 
for empirically accessible systems. Here, "empirically ac- 
cessible" is taken to indicate that the microscopic details 
of the system-bath interaction are available not merely 
in principle but in actual fact. (We might properly speak 
of an empirically accessible membrane, but this distinc- 
tion also will not bother us here.) The motivating (but 
by no means only) example is that of digital communi- 
cation networks (see, e.g., 0, 0), where a particular 
network interacting with a larger network (e.g., the In- 
ternet), constitutes the system-bath pair, and the data- 
exchange interactions are mediated by routers, switches, 
etc. that collectively act as the membrane between the 
system and bath. In such system administrator 

or researcher has access to everything that is happening 
at the membrane on a microscopic level: i.e., the indi- 
vidual packets being routed through the membrane are 
available for inspection by a passive demon. 

Before going further let us consider a more traditional 
example, that of a thermally conducting container filled 
with gaseous helium immersed in a bath of gaseous ar- 
gon. Here the membrane is the container itself, and a 
passive demon monitors it, keeping track of all the im- 
pacts of the atoms on the two surfaces of the container, 
and performing the relevant averages. It can be taken for 
granted that, when the substances are in thermal equi- 
librium, this demon is able to compute the temperature 
Jsj. The question therefore arises as to what makes the 
latter case any different in principle than the former, or 
some other one like a granular medium 

We claim that there is no essential difference. If, like 
our notional passive demon, we can obtain mesoscopic 



temporal and "spatial" averages as an observer at the 
membrane, and also keep records pointing to any infor- 
mation discarded in the averaging process, then we ought 
to be able to do thermodynamics from the bottom up, 
rather than from the top down. That is, an empirically 
accessible system affords us the opportunity to do ther- 
modynamics in order to describe the system. This is a 
stark departure from the usual cases, in which systems 
are empirically inaccessible and thermodynamics is used 
as a way of predicting the behavior of a system whose mi- 
crodynamical equations are known, but whose microstate 
is effectively unknowable. 

In the course of this paper we will demonstrate how 
empirical temperatures can be calculated for any finite 
system whose trajectory through state space is avail- 
able. As a corollary, we will obtain severe algebraic and 
topological restrictions for any empirical temperature de- 
fined as a function of occupation times. We will also 
describe an essentially universal mesoscopic data format 
for empirically accessible systems, enabling the explicit 
computation of temperatures. This in turn enables the 
entire thermodynamical/statistical- physical apparatus to 
be brought to bear on real-world problems. With that in 
mind, we will finally sketch an application. 



ISOTHERMS 

Most of this section derives from 0. We begin with 
some preliminaries. For the purposes of clarity and con- 
creteness, we will be fairly specific in some of the details, 
but much (if not most) of the discussion to follow ad- 
mits generalization/ variation on the minor points. We 
restrict the discussion to systems with a fixed, finite state 
space, where the states are labeled {1, 2, • • • ,n} for con- 
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TABLE I: Notations employed for the n— dimensional spaces 
used in this work. Both the vector notations and their coor- 
dinates are shown. 



Experimentalist 




Theorist 








{tl{Si), ■ ■ ■ ,t„{Si)} 




■■■ ,En-lis,),0{s,)} 



venience and have variable (but as yet undetermined) 
energies {Ei,E2, • ■ ■ , En}- The trajectory of the system 
through the state space evolves in continuous time start- 
ing at t = 0, and (again, for clarity and concreteness) 
the trajectory can be decomposed into disjoint epochs 
according to (for example) the following protocol. 

The (n -I- l)th epoch starts at i = ioo(so = 0) -I- ... -I- 
ioo(sra) smd ends as soon as every state has been visited 
at least once during the epoch: its duration— the Carlson 
depth of the ith epoch— is denoted ioo(sn)- 

In this way we will effect a mesoscopic moving time 
average over a fixed "length" scale (reflected in the cardi- 
nality of the state space) . During the ith epoch, the state 
occupation times are taken to be {ti(si), • • • ,i„(si)}, so 
that 



^tk{Si) = tc 



,{si 



(1) 



A second preliminary point deserves mention. Al- 
though the arithmetical structure of statistical physics 
allows the arbitrary numerical evolution of state energies, 
in a canonically distributed system we are constrained by 
the physics to interpret the underlying dynamics as si- 
multaneous changes to the energy levels and to the bath 
temperature. Another way of saying this is that temper- 
ature parametrizes the total bath energy. As a result, we 
are free to stipulate that the state energies always satisfy 



J2e, = o. 



(2) 



k=l 



This is just re-attributing homogeneous fluctuations in 
the state energies to temperature fluctuations. With this 
stipulation in mind, we introduce an empirical temper- 
ature parameter 0, and two sets of coordinate systems: 
experimentalist's and theorist's. See tabled 

It is natural to ask how to transform from one set of 
coordinates to the other. The transformation in one di- 
rection, from H— to space, is the conventional aim 
of statistical physics. The other direction is far less ex- 
plored, and negotiating it will allow us to bring the tools 
of thermodynamics and statistical physics to bear on a 
host of real- world problems. Without any constraints 
(e.g., the internal energy of the system is not fixed), a 
Boltzmann-Gibbs Ansatz for the state occupation prob- 



abilities is underdetermined, i.e., the assignment 



too (^i) 



(3) 



still does not uniquely specify a point in if— space but 
only a ray. 

Suppose that the state probabilities change slightly, 
say from one epoch to the next: we have the innocuous 
equation 



Ap["\si,Si+i) = Apf\s^,s^+i) 



(4) 



which embodies the essence of the Matched Invariants 
Principle (MIP): probability shifts agree in both coordi- 
nate systems. 

In what follows, we can dispense with the specification 
of particular epochs. The MIP implies that the total dif- 
ferential of each of the state probabilities can be written 
in two ways: 

I dt) = dpi; (V^pf ) I dH) = dpf (5) 

Now for each set of gradients, only rt — 1 are linearly in- 
dependent. In particular, the sets of gradients each span 
tangent spaces to spheres, because probabilities are con- 
stant on rays in both H— and t— space. Consequently, 
we can expand the differentials in terms of (say) the first 
n — 1 gradients: 



l|V«pf 'II 



drf ; dt 



n-l 

k=l 



\wk^\ 



drU6) 



Here we have introduced the differentials of new afhne 
parameters in the dr quantities. They are defined so that 
© holds. Because the first n — 1 gradients are linearly 
independent and span tangent spaces to spheres, these 
affine differentials are well-defined quantities. 

Combining jSJ and © immediately yields the thermo- 
dynamic evolution equation (TEE) 



dr 



H 



(7) 



where we have introduced the matrices 



(v*pi*^ I ) 



(H) 



HjH) 



IWf II 



jk 



\V"p 



(8) 



The TEE is a system of first-order ordinary differen- 
tial equations, and the standard existence and uniqueness 
results apply. The integral curves of the TEE in either 
space are circular arcs transverse to the foliation induced 
by the isotherms (i.e., the set of level hypersurfaces of 
constant temperature). This generic transversality con- 
dition implies, among other things, that the isotherms 
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are pinched at the origin in t— space. The Jacobian J 
simply relates the length parametrizations of these arcs 
under a change of variables, and we can use this fact to 
obtain isotherms numerically. 



EXAMPLE: A TWO-STATE SYSTEM 

In the case of a two-state system, the preceding consid- 
erations suggest working in terms of a single probability 
and a single energy, i.e., we have that 



Pi = P, P2 = 1 - p; Ei = E, E2 = -E 
and some routine algebra leads from here to 
1 i , i - p 



P ■■ 



1 



E 11 



tan( 



P 



(9) 



(10) 



The probability gradients are a bit more involved to 
compute: 



1 



to 



p 

1-p 



^"p-^{%'^) (11) 



The somewhat complicated appearance of the 
if— gradient is due to the constraint implied by the two 
probabilities summing to unity. In this example, the TEE 
is a scalar ODE, and we find that 



e 



lvy*)|| 



fF,2 



where 



IF, 2 



4 + log2(p-i - 1 



(12) 



(13) 



A geometrical argument (figures 1 and 2) allows for a 
semi-analytical determination of the isotherms in terms 
of a bijective uniform temperature map (UT map), abu- 
sively denoted 



e : {t : t = ^u'} ^ {H : H = Qu''} 

where we introduce the unit uniform vectors 
1 



u 



il; u 



H 



In 



ee 



{0,0,--- ,1}. 



(14) 



(15) 



The requirement for bijectivity of the UT map is just a 
codification of the zeroth law of thermodynamics. 

The semi-circles in iJ— space (see figure ^ map to 
the corresponding quarter-circles in space (figure EJ. 
Given a uniform temperature map, the concomitant 
isotherms in t— space can be identified numerically by 
integrating the TEE. 



00. 

'0. 



00, •/ 



FIG. 1: Schematic geometry of the transformation for a two- 
state system (higher dimensions are not qualitatively any dif- 
ferent). 




FIG. 2: Schematic of an isotherm in t— space. 



Briefly, suppose that we begin with the corresponding 
uniform points 

to - =Toti*; Ifo = e(ro)M^ = eo«^. (16) 

and jointly follow the corresponding integral curves ac- 
cording to 



dr* = To a((/)o,*); / dr^ = ©o 0o 



(17) 



The angles are defined in accordance with the figures, 
which are also important aids to the following arguments. 
We have that 0* = 0o cos((/)o,*), and combining the last 
few equations, we find that 



a('/'o,*) 



1 



To Jo 



00,. 



dr 



H 



dr 



H 



*-^0 



'0 "'O 



_ 60 j^" 

To Jo 
COS(' 



-jj d(j) 



dr 



too{(t>) 



Referring to the geometry of the H ^ t transforma- 
tion, we see that by integrating along circular arcs this 
could be solved numerically. However there is not much 
point in bothering: the differential equation will do just 
as well for numerical purposes. The point in this exercise 
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is twofold. First, it highlights the geometry governing the 
limits of integration for the TEE for an arbitrary number 
of states. Second, it details the need for, and the use of, 
a UT map in closing the transformation equations. 



PARAMETERIZATION OF THE THE 
ISOTHERMS 

From the preceding considerations it is clear that once 
a UT map is selected, the isotherms are obtained. Al- 
though the detailed geometry of the isotherms depends 
on the choice of UT map, some physical considerations 
allow us to immediately single out a particular class of 
candidate UT maps that all exhibit the same qualitative 
features. 

Let us return momentarily to the specific example of 
the container of helium immersed in argon. To a first ap- 
proximation, we can consider the atoms classically with 
conservative hard-body interactions. Imagine that in this 
idealization an active demon, through a homogeneous di- 
latation of the velocities, doubles the kinetic energies of 
all the argon atoms. The temperature (as the volume 
averaged kinetic energy of the atoms) that our passive 
demon computes is also doubled. 

Now let us assume that the passive demon has a very 
poor sense of inertia, and that our active demon does 
something else in addition to the foregoing: he replaces 
the passive demon's argon atoms with neon atoms of half 
the mass. The result of this trickery is, of course, that 
when the passive demon computes kinetic energies and 
temperature the original values are maintained. This 
suggests that if the same events happen more quickly, 
the temperature ought to be higher [6j. 

We can see as a result that the temperature ought to 
approach zero (by bijectivity, monotonically) as the Carl- 
son depth approaches infinity. Conversely, the tempera- 
ture ought to approach infinity as the Carlson depth ap- 
proaches zero. A UT map satisfying these requirements 
will be called topologically admissible (TA), as will the 
concomitant temperature function on t— space. Now we 
have a quite detailed characterization of any topologi- 
cally admissible temperature function. A caricature is 
provided by two examples below (Figures 3 and 4). 

It is perhaps natural to wonder if a particular TA UT 
map (or subclass of TA UT maps) might apply to every 
system. Although an answer might appear intractable, 
field-theoretical considerations (e.g., examination of the 
KMS condition 0) appear to suggest a UT map of the 
form Q = j^, where k is some constant. Publication 
of work in this direction is in progress by one of the au- 
thors and will be reported elsewhere, together with a 
straightforward algebraic derivation of the isotherms in 
terms of a UT map. 




FIG. 3: Temperature function for a two-state system with the 
uniform temperature map (j^)'' with 7 = 1. Both this figure 
and the next one were generated by an analytic formalism 
due to one of the authors that agrees with the numerical 
approach discussed here. 




FIG. 4: Temperature function for a two-state system with 
the uniform temperature map (j^)^ with 7 = 4. 



THE CLOSED QUEUE FORMAT 

Now that we have an empirical temperature, the is- 
sue arises of how to apply it to an empirically accessi- 
ble system. Towards that end we return once more to 
the metaphor of the passive demon monitoring abstract 
transactions across the membrane. In the case of our 
idealized gases, these transactions fall into several classes 
according to some discretization of the phase space. For 
instance, consider an argon atom with momentum p im- 
pinging upon the outside of the container at position x, 
where it donates some energy to the container and scat- 
ters. 

As a simplified energy transfer mechanism we might 
assume that the donated energy remains relatively local- 
ized as a phonon (it does not have time to diffuse much) 
before being absorbed by a helium atom impinging on 
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the inside. The passive demon can say that a transition 
of the form 



(p, x; Ar) 



phonon 



(p, x; He) 



(19) 



occurred. We might generahze this example somewhat 
so that both system and bath contain a mixture of aU 
the noble gases. 

Having discretized phase space near the membrane, the 
passive demon can associate triplets as above to discrete 
indices. The demon might choose to consider other ob- 
servables too. Then a transition of the form above is 
recast as 



. phonon '. 

3 » 3 



(20) 



In general, we require that all indices for the system and 
bath are common: say that there are B of them. The 
map from the tuples in H19|l to the indices in 1)20(1 defines 
the mesoscopic observables or bucket labels. 

Our mesoscopic data format is now described by fixing 
some positive integer b and considering the bucketspace 
or set of configurations of b indistinguishable balls in B 
buckets: 



as =b} 



(21) 



At the start of each epoch, the passive demon initializes 
a new mini-trajectory through bucketspace as follows: 
the trajectory begins at a fixed point in bucketspace. 
Every time a transition of the form above occurs, the 
passive demon moves a ball from bucket j to bucket j if 
such a move is possible: otherwise, he takes a ball from 
bucket j and replaces it. This sort of construction can 



trace its heritage back to the Ehrenfests |3| and Siegert 
0, among others. 

At the end of each epoch, the passive demon reports 
data like the transition rates and transition probabili- 
ties between buckets, the average state, etc. as well as 
the occupation probabilities Pa- Using these occupation 
probabilities, the temperature (and from it, other quanti- 
ties of interest) can be computed according to the scheme 
outlined above. More can be done: characterizing equi- 
libria, fiuctuations, and symmetries are all the subjects 
of ongoing work by the authors [loj . 

A sketch of how the preceding considerations can be 
fruitfully applied to real- world problems is in [lll | . 

The authors thank Debbie Huntsman for help in gen- 
erating figures, and James Luscombe for helpful conver- 
sations. 
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